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Our knowledge about the physical processes determining the activity of comets were mainly influenced by several extremely 
successful space missions (Giotto, Deep Space I, Stardust, Deep Impact and EPOXI), the predictions of theoretical models and 
the results of laboratory experiments. However, novel computer models should not be treated in isolation but should be based 
on experimental results and should be verified and calibrated by experimental work. Therefore, a new experimental setup was 
constructed to investigate the temperature dependent sublimation properties of hexagonal water ice and the gas diffusion through 
a dry dust layer covering the ice surface. We show that this experimental setup is capable to reproduce known gas production 
rates of pure hexagonal water ice. The reduction of the gas production rate due to an additional dust layer on top of the ice 
surface was measured and compared to the results of another experimental setup in which the gas diffusion through dust layers 
at room temperature was investigated. We found that the relative permeability of the dust layer is inversely proportional to its 
thickness, which is also predicted by theoretical models. However, the measured absolute weakening of the gas flow was smaller 
than predicted by models. This lack of correspondence between model and experiment may be caused by an ill-determination of the 
boundary condition in the theoretical models, which further demonstrates the necessity of laboratory investigations. Furthermore, 
the impedance of the dust layer to the ice evaporation was found to be similar to the impedance at room temperature, which means 
that the temperature profile of the dust layer is not influencing the reduction of the gas production. Finally, we present the results of 
an extended investigation of the sublimation coefficient, which is an important factor for the description of the sublimation rate of 
water ice and, thus, an important value for thermophysical modeling of icy bodies in the solar system. The achieved results of this 
laboratory investigations demonstrate that experimental works are essential for the understanding of the origin of cometary activity. 

Keywords: Comets, dust, Comets, dynamics, Comets, nucleus, Ices. 



(N 
> 
00 

in 



X 



1. Introduction 

Over the last decade, our understanding of the physical pro- 
cesses determining the activity of comets has deepened con- 
siderably. The main driving forces behind this progress were 
extremely successful space missions to comets: Deep Space I, 
Stardust, Deep Impact and, most recently, EPOXI. These mis- 
sions delivered many new and sometimes unexpected informa- 
tion, which urgently requires the development of new theoreti- 
cal models as well as new laboratory investigations. The con- 
struction of new computer models should not be treated in iso- 
lation but should be based on experimental results and should 
be verified and calibrated by experimental work. 

The most extended and challenging series of the comet sim- 
ulation experiments (known as KOSI) were performed in the 
1980s and 1990s using the German Space Agency's hardware 
facilities and the space simulation chamber in Cologne. The 
experiments dealt with the energy balance of insolated ice-dust 
mixture, the evolution of structure and composition of the sam- 
ple, ice sublimation and emission of dust particles. A review 
of the main results and their critical analysis can be found, for 
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example in Lammerzah l et al.| ( |1995| l and |Sears et al.| (fl999). 
Although the KOSI experiments provided new insights into the 
morphology and physical behavior of comet analogs, the com- 
plicated setup of the experiments significantly impeded a quan- 
titative analysis of the observations and did not cause a corre- 
sponding development of theoretical models. 

In order to avoid such complexity and to create a solid ex- 
perimental foundation for further theoretical investigations, we 
decided to carry out a set of small-scale experiments performed 
under conditions approximating those in the environment of 
cometary nuclei. The basic idea of these experiments is to use 
accurate methods for the investigation of the physical behavior 
of well-defined cometary analog materials. In this paper, we 
present the results of laboratory investigations on the sublima- 
tion of a solid hexagonal water ice and the reduction of subli- 
mation rate by a porous non-volatile dust layer covering the icy 
sample. 

Similar explorations have been executed in several labora- 
tories since the 1990s. Speaking about such experiments, one 
should first mention the extended series performed in the lab- 
oratories of Tel- Aviv University and the University of Graz. 
The Tel- Aviv group focused mainly on the study of amorphous 
water ice, its crystallization and the release of other volatiles 
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trapped by the amorphous water ice. In the recent publications 
( |Bar-Nun et"aLl|2008j F at ' E1 et al.||2009) , they presented fasci- 
nating results obtained for samples of gas-laden amorphous ice, 
which well reproduce the findings by the Deep Impact mission. 
The Graz group studied the energy and mass transport in both, 
solid and porous, crystalline water ice, including phenomena 
occurring when the icy sample was covered by a non-volatile 
layer. Their last experiments dealt with the exploration of the 
so-called solid-state greenhouse effect, which may play an im- 
portant role in the energy balance of icy surfaces in the solar 
system (Kaufmann et al. 2007| l. All these achievements are 



impressive. However, it is surprising that until now we have 
not obtained accurate systematic values of the effective subli- 
mation rate of water ice and the weakening of the sublimation 
beneath a porous dust layer. This paper is organized the fol- 
lowing way: in Sect. 2, we describe the experimental approach 
for the determination of the sublimation rate of hexagonal water 
ice and the impedance to the evaporation exerted by overlaying 
dust layers; Sect. 3 summarizes our experimental results for the 
sublimation of hexagonal water ice, the gas diffusion through 
dust layers, and the temperature dependence of the sublimation 
coefficient; finally, Sect. 4 indicates possible applications to icy 
bodies in the solar system. 



2. Experimental 




1: High-Vacuum Chamber 
2: Ice-Dust Sample 
3: Cold Plate 
4: Cooling System 



5: Temperature Sensors 
6: Pressure Gauge 
7: Chopper 
8: Aperture Plates 



9: Light Barrier 

10: Lock-In Amplifier 

11: Computer System 



Figure 1: Schematic diagram of the sublimation experiment. Sublimated wa- 
ter molecules from the sample's surface are detected using a pressure gauge 
together with a rotating chopper wheel. 



To investigate the temperature dependent sublimation prop- 
erties of hexagonal water ice and the gas diffusion inside a dry 
dust layer covering the ice surface, a new experiment consisting 
of an ice-dust sample under high vacuum condition and cryo- 
genic temperature was constructed (sublimation experiment) in 
our institute. The gas diffusion through identical dust layers 
was also measured with another experimental setup (diffusion 
experiment) in order to verify the obtained experimental results. 
This second experiment was performed at room temperature to 
avoid condensation and evaporation of gas molecules inside the 
dust layer. In this section, the technical setups, the experimental 
procedures and the sample preparation are described. 

2.1. Sublimation experiment 

The design of the sublimation experiment is shown in Fig. [T] 
The measurements with this setup were performed in a high- 
vacuum chamber (1) at pressures below 10~ 2 Pa. An ice-dust 



given in Sect. 2.4 The sublimated water molecules were col- 



sample (2, see Sect. 2.3 i was frozen onto a cold plate (3), which 
was connected to the cooling system (4). The entire system 
was cooled down to approximately ~ 110K with liquid nitro- 
gen (N2). Temperature sensors (5) inside the ice-dust sample 
and the cold plate enabled temperature measurements of the ice 
surface and the cooling system. With this experimental setup, 
the gas production rate of the ice-dust sample can be measured, 
using a pressure gauge (6) positioned 54 cm above the sample's 
surface and a rotating chopper wheel (7), which was used to 
periodically interrupt the flux of gas molecules. Therefore, the 
rotating chopper wheel was located between the sample and the 
sensor. A complete description of the gas flux measurement is 



limated into a narrow beam by two aperture plates (8) to elimi- 
nate the background gas flux. The signals of the pressure gauge 
and the light barrier (9), which was used to measure the rota- 
tional frequency of the chopper wheel, were input into a lock-in 
amplifier (10) to determine the sublimation rate of the ice-dust 



sample (see Sect. 2.4 1. Measurements can be performed auto- 
matically every two seconds by a computer system (11). The 
results of the data analysis are discussed in Sect. [3] 

In total, two different types of sublimation experiments were 
performed. For the calibration of the experimental setup, the 
gas production rate of hexagonal water ice was measured and 
compared with previous theoretical and experimental works 
(see Sect. 3.1 1. Afterwards, the ice-dust samples were used to 



investigate the gas diffusion inside dry dust layers covering the 
sample's surface and their influence on the sublimation proper- 
ties of water ice (see Sect |3.2| i. 

Before the start of a measurement sequence, the high-vacuum 
chamber was filled with nitrogen gas to avoid condensation 
of humidity onto the ice samples. Afterwards, liquid nitro- 
gen was led through the pipes of the cooling system in order 
to cool down the cold plate. At a temperature of ~ 255 K, 
a water droplet was used to freeze the sample onto the cold 
plate. Then, the high- vacuum chamber was evacuated from at- 
mospheric pressure to a pressure below 10~ 2 Pa. To exclude 
intermolecular collisions between the gas molecules, the mea- 
surements of the gas production rate and the temperature were 
performed below 1.14x 10~ 2 Pa, when the mean free path of the 
gas molecules, A, was larger than the distance between the sam- 
ple's surface and the pressure gauge. A detailed explanation of 
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the gas production rate measurements is given in Sect. 2.4 The 



thickness of the dust layer was the only parameter which was 
varied among the sublimation experiments. 

2.2. Diffusion experiment 

To investigate the gas diffusion through dust layers (see Sect. 
2.3 i with a different method, a second experiment was con- 
structed (diffusion experiment, see Fig. |2j. A dust layer of 
given properties (1) was positioned onto a filter paper (2), be- 
tween two pressure gauges (3) inside a vacuum chamber (4), 
which was evacuated to a pressure of ~ 10 _1 Pa by a turbo- 
molecular pump (5). A flow meter (6) enabled a precise control 
of the gas flux (7) through the dust layer. The dust layers in- 
vestigated with this experiment were not cooled to prevent con- 
densation and evaporation and to guarantee an investigation of 
the pure gas-permeability process. 




1: Dust Layer 

2: Filter Paper 

3: Pressure Gauges 

4: Vacuum Chamber 

5: Turbomolecular Pump 

6: Flow Meter 

7: Gas Flux 



Figure 2: Design of the diffusion experiment. The gas flux (dashed line) was 
forced to penetrate through the dust layer by the flow meter and the vacuum 
pump. 

For each experiment, four different gas fluxes (air) through 
the dust layer (pjf were adjusted with the flow meter: 0.56 x 
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1.31xl(T 3 Pam 3 s" 1 



and 



1 .69 x 10 Pa m s . These settings resulted in mean pressures 
ranging from 8.4 Pa to 105 Pa. After the gas flux adjustment, 
the pressure above and beneath the dust layer and the filter pa- 
per was measured and the permeability of the dust layer was 



computed (see Sect. 3.2 1. The permeability of the filter paper 
without dust layer was used for the normalization of the ob- 
tained results (see Sect. |3.2) . The results of these measurements 
are then compared with the obtained data from the sublimation 



experiment in Sect. 3.2 



2.3. Sample preparation 

In the sublimation experiments, cylindrical, hexagonal- 
water-ice samples with a diameter of 25 mm and a height of 
30 mm were used. The ice samples (see Fig. [3} were produced 




Figure 3: Pictures of two different cylindrical, hexagonal- water-ice samples 
inside the sample holder (1). For the temperature measurement, a thermocouple 
(2) was frozen into the ice close to the surface (a). Picture (b) shows a porous, 
dry dust layer composed of (54.7 ± 11,0) fan sized silica spheres, which was 
added on top of the ice sample (ice-dust sample). 



with distilled water that was frozen within a cylindrical sample 
holder at a temperature of ~ 255 K. To minimize the heat flow 
through the sample holder, a material with a low heat conductiv- 
ity (polyvinyl chloride, k = 0.15 WKr 1 rrT 1 ) and a small wall 
thickness (1 mm) in relation to its diameter was used. Due to 
the geometry of the sample holder, sublimation from the sam- 
ples's side walls was prevented. Both, the heat flow through the 
sample holder and the sublimation of molecules from the sam- 
ple's side walls, can have a considerable effect on the energy 
balance of the sample and, thus, on the surface temperature of 
the ice. For the temperature measurement of the icy surface, 
a thermocouple was frozen into the ice sample at a distance of 
< 3 mm from the surface. 

In eleven sublimation experiments (Sl-Sll, see Table |2| a 
dry dust layer composed of (54.7+ 1 1 .0) pm sized silica spheres, 
was added on top of the ice sample (ice-dust samples). To esti- 
mate the height of the dust layer, the mass of the dust layer was 
measured before and after the experiments. The size distribu- 
tion of these particles was measured with a microscope. Maxi- 
mum compaction of the dust resulted in a constant porosity of 
t// = (40.8 ± 0.3) % for the entire set of experiments. The poros- 
ity was determined by measuring the mass of water, which was 
filled into the voids of the dust. Furthermore, the compacted 
silica dust layer was also used in 17 diffusion experiments (Dl- 
D17, seeTable|3). 

To investigate the gas diffusion through a monolayer of dust, 
133 glass spheres with a monodisperse diameter of 2 mm, were 
positioned on the ice surface (SO) as well as on the filter paper 
(DO). This arrangement resulted in a porosity of i// = 43.4%, 
which is very close to the porosity of the dust layer. 

2.4. Measurement of the gas production rate 

The measurement of the gas production rate during the sub- 
limation experiments was conducted by using a Bayard-Alpert- 
type pressure gauge positioned above the sample's surface and a 
rotating chopper wheel, which was located between the sample 
and the sensor. Two apertures with a diameter of 25 mm were 
drilled into the chopper wheel, enabling the water molecules 
originating from the ice-dust sample to reach the pressure 
gauge. To increase the signal-to-noise ratio, the frequency of 
the rotating chopper wheel was chosen to less than the inverse 
response time of the pressure gauge. Conversely, high rota- 



Table 1: Results of different empirical measurements of the saturation pressure of water ice. The parameters ai and 02 refer to Eq.|3] 



Reference 



Fanale and Sal vail 



(1984) 



d\ [xl0 12 Pa] 



a 2 [K] 



Mauersberger and Krankowsky (2003 ) 
Mauersberger and Krankowsky ( 2003 ) 
This work 



T > 
T < 



169K 
169K 



3.56 
o 44 +0.09 

758.58 _22i 5 

-2 70 +6.73~ 
J J -0.38 



6141.667 
6132.9+ 1.8 

7044 ± 60 
6134.6± 17.0 



tional frequencies assure an improved signal stability. Thus, 
as a compromise, frequencies between 2.0 s _1 and 3.3 s _1 were 
used to meet both requirements. Using a rotating chopper wheel 
dissects the signal into two parts: in the first case, the rotating 
chopper wheel reflects all water molecules produced at the sur- 
face of the ice-dust sample; hence, only molecules of the sur- 
rounding gas are able to hit the sensor; in the second situation, 
molecules originating from the ice sample's surface and ambi- 
ent gas molecules can be detected by the pressure gauge. The 
pressure difference between both situations 6p was derived by 
an analysis of the modulated sensor signal with the lock-in am- 
plifier. The quantity 5p is equal to the pressure increase caused 
by the sublimating ice-dust sample. The rate of water molecules 
hitting the sensor j] s was calculated assuming elastic collisions 
between gas molecules and the sensor's surface A s , 



T]s(Td 



dp 



2m(v(Ti)) 



[»-]■ 



in 



Here, m is the mass and (v(r,)) the mean velocity of a water 
molecule, which depends on the temperature of the ice surface 
Ti. 

For comparison, the theoretical rate of water molecules hit- 
ting the sensor 77,, which was positioned at the distance r to the 
sample's surface, can be calculated by 



Tj t (Ti) = ff (r,)Z(r,)Q 



A, 



(2) 



Here, A,- describes the area of the sublimating ice surface and 
£1 = A s j r 2 ) is a factor that emerges from the geometry of the 
experimental setup. The ji r 2 term results from the integration 
over the cosine angular distribution of emitted gas molecules 
respective to the surface normal of the sample, according to 



Knudsen (1909a I. Due to the geometry of the sample holder, 



only molecules from the top of the cylindrical sample are able 
to sublimate. The sublimation rate Z(T{) of a pure and solid wa- 
ter ice surface is given by the classical Hertz -Knudsen formula 

|Knudsen| ( fl909b"l l 



Z(Ti) = p sub (Ti) 



[kgm 2 s 'J . 



2nkTi 

where k is Boltzmann's constant. The sublimation coefficient 

Z'(Ti) 



(3) 



Z(Td 



(4) 



describes the temperature dependent deviation of the measured 
sublimation rate Z'iTj) to the values predicted with the classi- 
cal Hertz -Knudsen formula (Kos sacki et al.| |T999) >. A detailed 
discussion of the temperature dependency of the sublimation 



coefficient is given in Sect. |3.3| According to the Clausius- 
Clapeyron equation and empirical measurem ents ([Fanale and 



Salvail 1984| Mauersberger and Krankowsky 2003 1, the satu 
ration pressure of water ice, p m b(Ti), can be formulated as 



PsubiJi) = ai e 



-a 2 /T, 



[Pa]. 



(5) 



Table [T] displays different obtained values for the coefficients 
a\ and a-i in Eq. [5] including the results of this work (see 
Sect. |3.1) , Alternative formulations of the temperature depen- 
dence of the saturation pressure of water are discussed by Buck 
( fl98T] l. 

3. Results 

3.1. Sublimation of hexagonal water ice 

For the calibration of the sublimation experiment, the gas 
production rate of hexagonal water ice rj s (Ti) was measured 
for different temperatures and compared with the theoretical 
gas production rate (see Eq. [2}. Fig. [4] shows the results of 
the measured gas production rates (open circles), as described 
in Sect. [2] The experiments were performed in a temperature 
range between 194.1 K and 151.7 K. For higher temperatures, 
the measurement of the gas production rate was not possible, 
due to the enhanced rate of evaporated water molecules and the 
associated pressure increase inside the high-vacuum chamber 
(see Sect. 2.1 1. Furthermore, the sublimation rate of hexagonal 
water ice above ~ 194 K cannot be explained with the classical 
Hertz-Knudsen formula, given by Eq. [3] (see sublimation coef- 
ficient, Sect. |3.3) , Thus, the calibration measurements started 
below this temperature. The lower detection threshold of the 
gas flux measurement was between 1 x 10 12 s _1 and 3 x 10 12 s _] 
and was given by the noise of the measuring system. The col- 
lected data were binned in intervals of 6 x 10~ 5 K -1 and the geo- 
metric means of the detected gas fluxes as well as the arithmetic 
means of the inverse temperatures were computed within these 
intervals (crosses). The statistical errors of these quantities are 
denoted by the error bars. For a comparison with the theoretical 
gas production rate of the ice sample (dashed curve, Eq. |2j, the 
binned data between 194.1 K and 163.9 K were logarithmically 
fitted to the function 



(6) 



(solid curve), with fitting parameters a\j and aip For these cal- 
culations, a sublimation coefficient of a,(r,) = 1 was used, 
because the investigated temperatures were below ~ 194K 



(see Sect. 3.3 1. The best agreement between the data and 
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Figure 4: Results of the gas production rate measurements of hexagonal water 
ice (open circles). The solid curve denotes the fit (see Eq. [6} to the geometric 
means of the binned data (crosses) between 194.1 Kand 163.9 K. For this calcu- 
lation a bin size of 6x 10~ 5 K _1 and a sublimation coefficient of ctj(T{) = 1 were 
used. A comparison with the theoretical gas production rate (dashed curve), 
given by Eq. [2] shows that the experimental setup is capable to reproduce the 
known gas production rates of hexagonal water ice over a range of temperatures, 
which were also measured by |Fanale and S alvail ( 1984 1 and Mauersberg er and] 
Krankowsky (2003 1. Note the deviation of the measured gas production rates 
to the theoretical curve below 163.9 K, which was caused by the noise of the 
measuring system. 



Eq. [6] was accomplished for ay = (3.23^^) x 10 12 Pa and 
a 2j = (6134.6 ± 17.0) K, which is in good agreement with the 
results from 



Krankowsky 



Fanale and Salvail ( 1984 1 and Mauersberger and 
(2003O" comparison of the obtained fit parame 



ters with this previous works can be found in Table[T] These cal- 
ibration measurements demonstrate that the experimental setup 
is capable to reproduce known gas production rates of hexago- 
nal water ice within the error tolerance. 

The major experimental error was caused by the uncertainty 
of the temperature sensor's position at the sample's surface. A 
displacement of 2 mm results in a maximal temperature error 
of 1 .6 K, which was comparable to the temperature uncertainty 
given by the manufacturer. This error can explain the scatter of 
measurements of the gas production rates. 

At temperatures below 163.9 K, the measured gas flux shows 
a systematic deviation from the theoretical curve, which was the 
result of the noise of the measuring system. Thus, we consider 
a temperature of 163.9 K as the lowest temperature accessible 
to our experimental setup. 

3.2. Gas diffusion through dust layers 

In total, the gas production rate of twelve ice-dust samples 



rjid(Ti, h) (see Sect. 2.3 i were investigated with the sublimation 
experiment (SO - Sll). The individual measurements are sum- 
marized in Table [2] The height of the dust layer h is given in 
units of the mean particle diameter. 

From the deviation of the measured gas production rates to 
the values for pure hexagonal water ice without a dust layer (see 
Sect. 3.1 1, the normalized gas production rates of the ice-dust 
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Figure 5: Example of the determination of the normalized gas production rate 
of an ice-dust sample in the sublimation experiment. The open circles display 
the measured gas production rates of experiment S9, for which the height of the 
dust layer was (56.56 ± 1.12) particle diameters. Fitting the geometric means 
of the data (crosses) below 194.2 K with Eq.|7J(solid curve) and comparing this 
result with the gas production rate of hexagonal water ice without an additional 
dust layer (see Eq. [6j gives the normalized gas production rate of the ice- 
dust sample. For comparison, the theoretical gas production rate of a pure 
ice sample is visualized by the dashed curve (see Eq. |2J. The best match 
to the data was achieved for B^y = (0.38^9|S) X 10 I2 Pa, which results in a 
normalized gas production rate of 0.118^ JJjJ'. An extrapolation of the fit to 
higher temperatures (dash dotted curve) demonstrates the strong variation of 
the sublimation coefficient above ~ 194 K (dotted line). For this data analysis a 
bin size of 6 X 10~ 5 K - ' was used. 



samples, ?7,v/(r,, /i)/?/,(7',), were calculated. For an explanation 
of this derivations, the results of the experiment S9 are shown 
in Fig. |5] 

In this experiment a dust layer with a height of (56.56 +1.12) 
particle diameters was used. The measured gas production rates 
(open circles) were binned in intervals of 6 x 10~ 5 K _1 and the 
geometric means of the gas production rate as well as the arith- 
metic means of the inverse temperatures within these intervals 
were calculated (crosses). The error bars denote the statistical 
errors of the derived values. To compare the gas production 
rates of the ice-dust sample with the gas flux of a pure (i.e. 
dust-free) hexagonal ice samples (dashed curve, Eq. |2j, the fol- 
lowing function was used to fit the data, 

rjid^h) = a(Td a ud (h) e-"^ /T ' Q. - [s^l (7) 

(solid curve in Fig. |5J. An additional dust layer on top of the 
ice surface decreases the gas production rate of the ice sam- 
ple. However, this reduction should not depend on the tem- 
perature of the sample, if we assume that the gas molecules 
are only scattered by the individual particles of the dust layer. 
If the gas molecules were absorbed on the surface of the dust 
particles, evaporation of the trapped gas molecules occurs over 
time. In this case, the rate of evaporation of the condensed wa- 
ter molecules inside the dust layer is bigger than the rate of 
water molecules supplied by the ice surface, if the temperature 
of the dust layer is higher than the temperature of the ice sur- 
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Table 2: Summary of the experimental results achieved with the sublimation experiment. The normalized gas production rates of the ice-dust samples 
Vid(Tj,h)/r)i(Ti) were calculated with Eq. [8] In the experiment SO, a monolayer of glass spheres with a diameter of 2 mm was positioned on the ice surface 
(see Sect. [23). 



Experiment 


h (Particle Diameter) 


aud(h) [xl0 12 Pa] 


nid{Ti,h)lt]i(Ti) 


S0' f 


1.00 ±0.00 


not +0.04 
Z -°-> -0.04 


0.882 


SI 


10.62 ±2. 12 


1 30 +oos 

1JU -0.07 


0.402 


S2 


17.10 + 0.11 


1 01 +008 

lul -0.08 


n t i q +0.024 


S3 


18.22 ±0.34 


96 +0 07 

-0.07 


n 907 +0.020 

v-^yi _ 00i9 


S4 


24.59 ±0.89 


0.76 


235 +0 015 

-0.014 


S5 


26.04 ± 1.68 


0.68 tlj 


211 +0031 

-0.027 


S6 


30.51 ±0.11 


63 +0 05 

uu -' -0.05 


195 +om5 

U - ,7J -0.013 


S7 


36.55 ±1.90 


56 +004 

UJU -0.04 


0.173^ 


S8 


47.62 ±0.89 


0.35 Tol 


0-108 ™> 


S9 


56.56±1.12 


0.38 


0.118 


S10 


60.13 ±1.34 


27 +0 05 

-0.04 


0.084^3" 


Sll 


67.85 ± 1.12 


25 +0 03 

"" iJ -0.03 


0.077 3$ 



f: Investigation of a monolayer of glass spheres instead of the dust layer. 



Table 3: Results of the diffusion experiments. The normalized rates of molecules penetrating through the dust layer r]jf(h)/T]f were computed with Eq. |18| using 
the measured permeabilities of the dust layers &j(h). In experiment DO, a monolayer of glass spheres was investigated instead of a dust layer. The experiments D2, 
D10 and D17 were extended to higher pressures to demonstrate that the measurements were not influenced by the increased pressure (see Fig. [5). 



Experiment 


h (Particle Diameter) 


@ d (h) [xlO- 6 m 3 s-»] 


ridf(h)/T] f 


DO 1 


1.00 ±0.00 


1514.59 ±240.71 


0.770 ±0.085 


Dl 


12.72 ± 1.12 


86.67 ±5.56 


0.342 ±0.021 


D2* 


13.50+ 1.12 


63.66 ±1.59 


0.393 ±0.006 


D3 


14.06+ 1.12 


41.01 ±8.01 


0.236 ±0.021 


D4 


19.19 + 1.12 


47.64 ± 10.75 


0.263 ±0.025 


D5 


20.53 ± 1.12 


58.00 ±3.24 


0.258 ±0.021 


D6 


24.55 ± 1.12 


26.92 ±4. 11 


0.215 ±0.021 


D7 


27.90± 1.12 


40.52 ±4.47 


0.1 95 ±0.007 


D8 


28.79 ± 1.12 


42.59 ± 1.87 


0.204 ±0.021 


D9 


32.36+ 1.12 


22.58 ±2.56 


0.187 ±0.015 


D10* 


34.70+ 1.12 


21.22 ±0.65 


0.177 ±0.004 


Dll 


35.04± 1.12 


34.78 ± 1.24 


0.173 ±0.016 


D12 


41.29 + 1.12 


21.99 ± 1.51 


0.1 16 ±0.008 


D13 


44.19 + 1.12 


18.88 ±1.29 


0.131 ±0.011 


D14 


50.43 ± 1.12 


19.29 ± 1.92 


0.133 ±0.005 


D15 


54.45 ± 1.12 


10.50 ±1.05 


0.077 ±0.001 


D16 


57.58 + 1.12 


13.55 ±1.06 


0.098 ±0.004 


D17* 


72.75 ± 1.12 


10.42 ±0.51 


0.096 ±0.004 



f: Investigation of a monolayer of glass spheres instead of the dust layer. 
Measurements extended to higher pressures. 
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face. Thus, the dust layer loses the water ice after some period 
of time. Even in this case, the influence of the dust layer on 
the gas production rate of the ice sample is still not affected by 
the temperature of the ice surface (a detailed discussion of the 
thermal influence of the dust layer on the gas production rate 
is given at the end of this Section). Thus, the temperature of 
the dust layer and, therewith, the temperature of the ice surface 
will not affect the influence of the dust layer on the gas produc- 
tion rate. Thus, we assume that the dust layer on top of the ice 
sample does not influences the slope of the function in Eq. [JJ 
Therefore, the parameter a2,irf> which determines the slope of 
the fit was fixed to fit the data, i.e. a^M — a 2j ( see Sect. [XT}. 
Thus, the coefficient a\ ^{h) was the only free parameter, which 
was varied to fit the binned gas production rates with Eq. [7J 
From the results of the fits, the normalized gas production rates 
(rates of molecules penetrating through the dust layer) were cal- 
culated, 



rndjj, h) _ aijd(h) 



(8) 



For the S9 experiment, the best match to the binned data was 
found for a lJd = (0.38+°^) x 10 12 Pa. This results in a normal- 
ized gas production rate of 0.1 18^'^ ■ 

The fit to the data could only be realized below a critical 
temperature of 194.2 K (solid line), due to a change in the 
sublimation coefficient at higher temperatures (see Sect 



3.3 1. 



To demonstrate the decrease of the sublimation coefficient at 
higher temperatures, the fit to the binned data was extrapolated 
(dash dotted line). 

The obtained fit parameters and the derived normalized gas 
production rates for the different sublimation experiments are 
presented in Table [2] For the SO experiment, a monolayer of 
133 glass spheres was added on top of the ice surface (see Sect. 
|2.3| >. The result of this experiment is visualized in Fig. [6] Note 
the deviation of the measured gas production rates to the ex- 
trapolated fit (dash dotted curve) above ~ 194 K, due to the 
decrease of the sublimation coefficient at higher temperatures 
(see Sect. [33}. 

Further investigations of the gas diffusion through dust lay- 
ers at room temperature were carried out using the diffusion 
experiment. In 18 experiments (DO - D17), the permeability of 
the dust layer and the normalized rate of molecules penetrating 
through the dust layer were measured at room temperature. The 
results of these experiments are summarized in Table [3] 

The Knudsen number, defined by 



Kn = 



*eff 



(9) 



with the mean free path of the gas molecules, A, and the effec- 
tive pore size of the material, d e /f, is characteristic for the flow 
regime. Free molecular flow is expected for Kn > 1. Due to 
the low gas pressures in our experiments, the Knudsen num- 
bers were always Kn > 3.3 so that collisions of gas molecules 
with the particular medium dominate over intermolecular colli- 
sions, and the mass is transferred due to the pressure gradient 
according to Knudsen's law, where gas molecules migrate in- 
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Figure 6: Results of the experiment SO, in which 133 monodisperse glass 
spheres were positioned on the ice sample to investigate the normalized gas 
flux through a monolayer of dust. The achieved data (circles) were binned in 
intervals of 6 X 10~ 5 K~' and the geometric means of the gas fluxes and the 
arithmetic means of the inverse temperatures were derived (crosses). The sta- 
tistical errors of thesesquantities are visualized by the error bars. The best fit 
to the data using Eq. 7 was achieved for ai t a = (2.85+jj^) X IO 12 Pa (solid 
curve). This results in a normalized gas flux of 0.882+jj 'jSI for the monolayer 
dust. Note the deviation of the detected gas fluxes to the extrapolated fit (dash 
dotted curve) at temperatures above 1 94.2 K. 



dependently of each other. The diffusive flux of an ideal gas is 
described by Fick's first law, 



dC r o 1 1 

Jmol = D Kn — |mol m s J 



dx 



(10) 



relating the diffusion flux j mo [ to the gradient of the molecular 
concentration C [molm 3 J. D^ n [m 2 s _1 ] is the Knudsen dif- 
fusion coefficient or diffusivity. With the ideal gas law, Fick's 
first law can be formulated in terms of the pressure gradient, 



Jmol 



D Kn dP 
R T g dx 



(ID 



where R is the ideal gas constant and T g the temperature of the 
gas. This equation can be rewritten for the flux j measured 
in numbers of molecules flowing through a unit area per time 
interval, 



Dk„ dP 
kT g dx 



[m-T 



(12) 



For the steady gas flow through a random porous medium, the 
effective Knudsen diffusivity can be introduced by 



l <A 

D Kn - -~d eff ~ 



8/tr„ 



(13) 



where if/ is the transport porosity of the medium and t is the 
tortuosity factor. The medium used in the diffusion experiment 
has only open pores. Thus, the transport porosity equals the to- 
tal porosity. The tortuosity factor is defined as the square of the 
ratio of the effective path (or actual length L e ) to the thickness 
of the porous layer h, or as the square of the ratio of the average 
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free path traveled by the species (/) to the average displacement 
in the direction of diffusion (x) ( |Boudreau| [1996 ) 

= % = Hi 
T h 1 (x 2 ) ' 



(14) 



From the measured pressure difference Ap and the adjusted gas 
flux through the dust layer and the filter paper <pdj, the perme- 
ability of the dust layer ®d(h) was determined by 



Ap 

4>df 



1 



[mV 1 ]. 



(15) 



Here, 0/ is the permeability of the filter paper, which was es- 
timated in several calibration experiments. This inter-relation 
arises from the fact that the permeability is given by the inverse 
resistivity to the gas flow through the material. This resistiv- 
ity can be derived from the quotient of the pressure difference 
(potential difference) and the gas flux (current). Thus, the dust 
layer acts like a resistor to the gas flux. In the case of the sub- 
limation experiment, the ice surface of the sample possesses a 
certain internal resistance, which is given by the inverse sub- 
limation probability of the water molecules. Furthermore, the 
internal resistance in case of the diffusion experiment is caused 
by the filter paper. Therefore, the resistivity of the dust layer 
together with the internal resistivity can be described with two 
resistors in series-connection (resistivity model). 

Fig. [7] shows the derived permeabilities as a function of the 
height of the dust layer (crosses). From the fit to these data 
(solid curve), the following correlation between the permeabil- 
ity and the height of the dust layer was found, 



®d(h) 



©o 
h 



(16) 



where O = (4.96 ± 0.98) x Kr^nrs -1 is the permeability 
coefficient of the dust layer. The result of the DO experiment 
was not included into the fit, to achieve a better match to the 
derived permeabilities in the interval between 12.72 and 72.75 
particle diameters. The dashed curve demonstrates the extrapo- 
lation of the applied fit. The major error in these measurements 
are caused by the uncertainty of the pressure signal and the gas 
flux adjustment. 

A comparison between Eqs. [T2] [15] and [TS] shows that the 
effective Knudsen diffusivity is connected with the permeability 
coefficient by the evident ratio 



D 



Kn 



©o 

A d 



(17) 



where A,/ is the cross section of the sample. Because the gas 
temperature and porosity of the sample are known, one can 
evaluate the ratio of the effective pore size to the tortuosity fac- 
tor. Further evaluation requires the use of a computer model de- 
scribing molecular diffusion in close packed porous media. The 
construction of such a model is in progress and corresponding 
quantitative analysis of the results will be presented in a forth- 
coming paper. 

To compare the results of the diffusion experiment with the 
results of the sublimation experiment, the normalized rates of 
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Figure 7: Derived permeabilities for different heights of the dust layer, calcu- 
lated with Eq. |15| (crosses) and fitted with Eq. |16| (solid curve). The result of 
the DO experiment was not included into the fit, to achieve a better match to the 
data between 12.72 and 72.75 particle diameters. From the fit, the permeability 
coefficient O = (4.96 ± 0.98) X 10~ 8 m 4 s _1 was derived. 



molecules penetrating through the dust layer and the filter paper 
rjdfQi)Mf were calculated using the resistivity model 



J]df(h) 
Vf 



®d(h) 



® f + ® d {h) 



(18) 



Fig. [HJpresents the comparison between the normalized gas pro- 
duction rates of the ice-dust samples (squares, Eq. [8} and the 
resulting normalized rates of molecules penetrating through the 
(warm) dust layer and the filter paper (crosses, Eq. 18 1. The er- 
ror bars indicate the statistical errors, which were derived from 
the errors of the fits (sublimation experiment) and from the stan- 
dard deviations to the arithmetic means (diffusion experiment). 
The decrease of the normalized diffusion rates for increasing 
heights of the dust layer were fitted using the resistivity model 



T]id,df(h) 



l 



1 +b- 1 h 



(19) 



Here, b is the height in particle diameters at which 50% of 
the initial gas flux penetrates through the dust layer. The best 
match to the normalized diffusion rates was found for b,d = 
(7.31 + 0.20) particle diameters in case of the ice-dust samples 
(solid curve) and bdf - (6.54 ± 0.38) particle diameters for the 
dust layers at room temperature. Both data sets show the same 
functional behavior with respect to the height of the dust layers. 
This accordance can have two different causes: 

1. Due to the low heat conductivity of the loose material 
(Krause et al. submitte d to Icarus) , the temperature de- 
crease rate inside the dust layer was much smaller than the 
cooling rate of the ice. Therefore, the temperature of the 
bulk of the dust layer remained high enough so that no con- 
densation and recondensation of water molecules occurred 
inside the dust layer. A comparison of the he at conductiv- 
ity of hexagonal water ice (3.26 WirT 1 K _1 ; Slack 1980 1 
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Figure 8: Derived normalized rates of molecules penetrating through the dust 
layer, measured with the sublimation experiment (squares) and the diffusion 
experiment (crosses) as a function of the height of the dust layer. The data were 
fitted with Eq. |19| (resistivity model). The best match to the data was achieved 
with bij = (7.3 1 ± 0.20) particle diameters for the sublimation experiment (solid 
curve) and bdf = (6.54 ± 0.38) particle diameters for the diffusion experiment 
(dashed curve). The height of the dust layer is given in units of the monomer 
diameter. Mind that two different monomer sizes were used in the experiments, 
i. e. 2 mm (monolayer experiment) and (54. 7± 1 1 .0) /im (any other experiment). 



with the heat conductivity of relat ively densely packed sil- 
ica particles (~ 10~ 2 W m" 1 K" 1 ; |Krause et al.| |submitted| 



to Icarus) demonstrates that the temperature gradient in- 
side the dust layer should be much higher than the temper- 
ature gradient within the ice sample. Thus, only the first 
few layers of the dust particles above the ice surface can 
remain cold enough for condensation of water molecules. 

If we assume that part of the dust layer was principally 
cold enough for ice condensation to take place, the evap- 
oration rate of water molecules from inside the dust layer 
must be, however, higher than the evaporation rate of the 
icy surface, due to the higher temperature of the dust 
layer. Thus, the dust layer cannot quantitatively store wa- 
ter molecules in form of icy condensates on the dust parti- 
cles' surfaces. 



In our experiments we measured the surface temperature of 
the ice surface, therefore, the thermal influence of the presence 
or absence of the dust layer is unimportant for the interpretation 
of the gas production rate of the ice surface. The temperature 
of the ice surface is influenced by back-scattered molecules in- 
side the dust layer and by the transported heat through the dust 
layer and through the ice sample, but the gas production rate is 
determined by the actual surface temperature, which was con- 
tinuously measured. However, in future experiments we will 
concentrate on the heat transport through the dust layers by in- 
solation. The influence of the grain size and the porosity on the 
heat transport through dust layers will also be studied. 

The match between the two experiments also implies that the 
diverse pressure regimes of the two different experimental en- 
vironments have not affected the measurements. In both exper- 



Figure 9: Extended permeability measurements for lower Knudsen numbers 
performed with experiments D2 (pluses), D10 (triangles) and D17 (squares). 
The means of the individual measurements are denoted by the dotted, dashed 
and dash dotted lines. 



iments, the mean free path of the gas molecules was larger than 
the mean pore size of the material, which was of the same order 
as the mean particle size. A calculation of the Knudsen num- 
bers Kn demonstrates that the experiments were performed in 
the free molecular flow regime, i.e. Kn = 9877.5 (sublimation 
experiment) and Kn = 3.3 (diffusion experiment), respectively. 
Thus, the obtained results should not be influenced by the differ- 
ent pressure regimes as long as the Knudsen numbers guarantee 
a free molecular flow through the dust layer. 

Further investigations with the diffusion experiment at higher 
pressures (D2, D10, D17) have shown that the permeability re- 
mained constant even for smaller Knudsen numbers (see Fig. 
[5J. For the calculation of the Knudsen numbers, the mean par- 
ticle diameter was used to determine the characteristic length 
of the system. Furthermore, the presence of a laminar flow 
through the dust during the diffusion experiments was con- 
firmed by an analysis of the Reynolds number, which was in 
the range between 0.8 and 9.8. The gas flow through the dust 
layer on top of the ice sample was obviously not influenced by 
turbulence. 

3.3. Sublimation coefficient 

The sublimation coefficient (see Eq. Q is introduced to de- 
scribe the temperature dependent deviation of the measured 
sublimation rate to the expected theoretical value of hexagonal 
water ice, given by the classical Hertz-Knudsen formula (see 
Eq. [3]and Fig. From the standpoint of statistical physics, the 
sublimation of water ice should be considered within the gen- 
eral theory describing adsorption and desorption processes. In 
the case when the adsorbate and adsorbent are identical and the 
vapor pressure is negligible in comparison with the saturation 
pressure determined from the Langmuir adsorption isotherm, 
the equation describing the sublimation of the ice crystal is re- 
duced to Eq. [3] which is widely used in planetary physics. 
It is common knowledge that this formula is trivially derived 
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from the kinetic theory of gas using just one assumption about 



the equilibrium distribution of molecular velocity ( Skorov and 



Rickman |1995[ ). This means that all information about the mi- 
crophysical processes determining the actual sublimation rate 
as well as the growth rate of the ice crystal is hidden (masked) 
in the so-called sublimation coefficient, which obviously has a 
complex physical background. A detailed microphysical inves- 
tigation of this coefficient is indeed a sophisticated theoretical 
task, which is far beyond the scope of this work. It should be 
only mentioned that the rate of sublimation is influenced by 
the shape of the crystal, the structure of its surface, the pres- 
ence of impurities and other physical factors. For example, 
the theoretical study of the sublimation rate for the case of an 
ideal homogeneous surface was presented by Shul'Man ( 1972] > 
where he treated the sublimation as a variant of a more general 
process of adsorption-desorption using the method of chemical 
potentials. A more complete consideration of the evaporation 



problem was carried out in the work performed by|Knake and 



Stranskii (1959) in which the authors took a set of parameters 



into account characterizing the heterogeneity of the crystal, i.e. 
a variable coordination number and, therefore, variable activa- 
tion energy of the molecules, the geometry of the crystal lattice, 
and the distance between the defects. It should be emphasized 
that it is the complexity of the problem, which makes the the- 
oretical formula unsuitable for practical use, because many of 
these microscopic characteristics are not well defined in labo- 
ratory experiments, and in addition, they can vary over time. 

This is why new direct experimental measurements 
of crystalline-ice sublimation/growth under particular non- 
equilibrium conditions (mainly given by the temperature and 
pressure range) are so important for planetary applications. The 
assumption that sublimation is the exact reverse of a growth 
process allows us to use experimental data obtained during the 
investigation of the latter (Pruppacher and Klett 19971. The 
corresponding set of references can be found in the work of 
|Kossacki efaX] ( |1999| ) and |Nelson| ( |1998| ). It should be noted 
that the deviation of the sublimation coefficient from unity were 
observed by Kramers and Stemerding (1951), Isono and Iwai 



( 1969 ), Beckmann and Lacmann ( 1982|, and later by Kossacki 



|et al.| ( fl999| ) with quite different experimental setups. At the 
same time, this effect was clearly observed at about the same 
temperatures, i.e. above 190K. At lower temperatures, there 
is a good agreement between the theoretical estimate from the 
Hertz-Knudsen formula (see Eq. [3} and the achieved experi- 
mental results. Hereafter, we present the results of our system- 
atic evaluation of the sublimation coefficient carried out under 
experimental conditions reasonable for the planetary science 
applications. 

The dashed curve in Fig. 



10 visualizes the temperature de- 



pendence of the sublimation coefficient, derived by Kossacki 



et al. ( 1999 1. This retrieved functional behavior was only based 



on three different experimental wo rks ([Beckmann and Lacmann 



([1982 squares), |Kossacki et al.| ( |1999 cross), and Davy and 



Samarjai ( 197 1 [ triangle)), performed in three limited temper- 
ature intervals. For further comparison, the results of previ- 
ous works are also presented in Fig. 
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Figure 10: Temperature dependence of the sublimation coefficient of hexagonal 
water ice, calculated from the gas production rates measured with the sublima- 
tion experiment (experiments: SO - S 1 1). The retrieved data (open circles) were 
binned in intervals of 10~ 5 KT 1 , and the geometric means of the derived subli- 
mation coefficients as well as the arithmetic means of the inverse temperatures 
were computed (pluses). According to these results, a new function for the tem- 
perature dependence of the sublimation coefficient is proposed (solid curve, see 
Eq. |20) . A good match to the data was retrieved for Ci = (1 1.58± 0.24)x 10 3 K 
and c 2 = (4.99 ± 0.20) X 10"" 3 K -1 . For comparison, the results of the fol- 
lowing works are also visualized: [Kramers and Stemerding \ 1951 1 diamonds) 
Isono and Iwai ( 1969 asterisks, red dash dotted curve), [Beckmann and Lac- 
mann ( 1982 squares), Davy and Samarjai ( 1971 triangle) and Kossacki et al 
1 1999 cross, dashed curve). The red dash dotted lines denote the expected 
temperature range of the sublimation coefficient measured by Isono a nd Iwai| 
(T969[ . 



However, the disagreement between these different experimen- 
tal works is evident. Thus, a more detailed investigation of 
the temperature dependence of the sublimation coefficient in a 
broader temperature range and under well defined experimental 
conditions was conducted, using the sublimation experiment. 

Using the results of the experiments SO - Sll, the sublima- 
tion coefficient was derived from the deviation of the measured 
gas production rates to the applied fits (see Eq. [7] sublimation 
experiment) and the related extrapolations (see Fig. B). Fig. 10 
presents the calculated sublimation coefficients (open circles). 
The data were binned in intervals of 10~ 5 K _I and the geomet- 
ric means of the derived sublimation coefficients as well as the 
arithmetic means of the inverse temperatures were computed 
(pluses). This calculations were only conducted in intervals, 
where three or more individual measurements were performed. 
Due to the retrieved results, a new formulation for the tempera- 
ture dependence of the sublimation coefficient is proposed, 



a{T t ) = 



0.854 



1 + exp[ - Cl (rr 1 - c 2 )] 



0.146 



(20) 



(solid curve), with ci = (11.58 + 0.24) x 10 3 K and c 2 = 
(4.99±0.20)xl0~ 3 K _1 . This new function is in good agreement 
with the results from Kossacki et al. (1999) for temperatures 
above ~ 227 K and below ~ 194 K. Between these two temper- 
atures, the new formulation deviates from the functional behav- 



merding ( 195 1 diamonds) and Isono and Iwai ( 1969 asterisks). ior published by Kossacki et al. ( 1999 ) (dashed line). However. 
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the new proposed temperature dependence of the sublimation 
coefficient is based on 4, 048 data points, which were achieved 
by eleven individual measurements. 

The different formulations of the temperature dependencies 
of the sublimation coefficient have a huge influence on the en- 
ergy balance of sublimating ice surfaces and thus on the thermal 
modeling of icy bodies in the solar system. 

4. Conclusions and implications 

The overt lack of experimental data on ice sublimation un- 
der non-equilibrium conditions (low pressure) and on transport 
properties (i.e. on the absolute permeability) of porous media 
results in a lack of understanding of cometary activity. There- 
fore, we focused on two principal questions in this paper: (1) 
how can sublimation be kept alive beneath the covering porous 
dust layer and (2) how does the sublimation rate depend on the 
temperature of the sublimating ice? 

We found that the normalized rate of penetrating molecules 
through a dust layer is inversely proportional to its thickness. 
This result can be explained with the introduced resistivity 



model (see Sect. 3.2 1, in which the dust layer is treated as a 
resistor to the gas flux. Using this model, we are able to calcu- 
late the thickness at which 50 % of the initial gas flux penetrates 
through the dust layer: (7.31 + 0.20) particle diameters (subli- 
mation experiment) and (6.54 ± 0.38) particle diameters (dif- 
fusion experiment). This accordance among the two different 
experiments may have two different causes. Either the temper- 
ature decrease rate inside the dust layer was much smaller than 
the cooling rate of the ice, due to the low heat conductivity of 
the loose material (Kra use et al.| [submitted to Icarus| l. There- 
fore, the temperature of the dust layer probably remained high 
enough so that no condensation and recondensation of water 
molecules have occurred inside the dust layer. Or, the tempera- 
ture of part of the dust layer was principally cold enough for ice 
condensation to take place. In this case, the evaporation rate of 
water molecules from inside the dust layer must be, however, 
higher than the evaporation rate of the icy surface, due to the 
higher temperature of the dust layer. Thus, the dust layer cannot 
quantitatively store water molecules in form of icy condensates 
on the dust particles' surfaces. 

The second objective of this work was the systematic evalua- 
tion of the so-called sublimation coefficient, defined as the ratio 
between experimental and maximal sublimation rate (Eq. ffl. 
Therefore, we demonstrated that the experimental setup is ca- 
pable to reproduce the known gas production rates of hexagonal 
water ice. From the deviation to the theory, the sublimation co- 
efficient was calculated for temperatures ranging from 223 .4 K 
to 157.4K. This investigation revealed that the sublimation co- 
efficient is unity for temperatures below ~ 194 K. Between 
~ 194K and ~ 227 K, the sublimation coefficient decreases 
by approximately one order of magnitude. This behavior has a 
tremendous influence on the energy balance of sublimating ice 
surfaces and, therefore, on the thermal modeling of icy bodies 
in the solar system. 

Until recently, the questions which are the subject of this 
work had only an indirect relation to the physics of comets, be- 



cause the existence of a non-volatile porous crust on the surface 
of the cometary nucleus as well as the sublimation of volatiles 
through this crust remained unproven hypotheses and theoreti- 
cal assumptions. However, the situation has changed dramati- 
cally over the past few years, when successful space missions 
(Deep Space I, Stardust, Deep Impact, EPOXI) delivered a lot 
of new facts relating to the physics of cometary nuclei. 

Today, we possess high-resolution images of five nuclei of 
periodic comets. These images, together with spectral and 
infrared observations, have given us solid evidence that all 
cometary nuclei observed till now are largely covered by a 
porous, dark, non-volatile crust. Especially exciting results 
were obtained by the Deep Impact mission to comet Tempel 1 
and later to comet Hartley 2. The thermal inertia of the surface 



of comet Tempel I was found to be very low (|Groussin et al. 



2007). This result indicates that a high-porosity material forms 
the surface layer (Krause et al. submitted to Icarus i. At the 
same time, only traces of water ice were found on the surface 



(Sunshine et al. 2007 1. These findings, together with the mea- 
sured gas activity ( Schleicher et al. 2006) can be considered as 
a strong argument in support of ice sublimation beneath a thin 
dust layer. The latest images of comet Hartley 2, taken during 
the flyby of the Deep Impact spacecraft (http://epoxi.umd.edu), 
give evidence of the high activity of the nucleus, with numerous 
jets visible on the surface. The jets are still alive on the night 
side of the nucleus. As in the case of other comets, the surface 
of Hartley 2 is dark, an indication that it is covered by dust. 
Thus, one can expect that the transport of sublimation products 
through a porous non-volatile dust layer is an important pro- 
cess on comet Hartley 2. All these data make the investigation 
of mass transport through porous media an interesting and im- 
portant endeavor. 

Recently, we discussed in detail the theoretical aspect of the 
first question in application to the cometary physics ( |Skorov| 
et al. submitted t o Icarus| . In that paper, we reviewed the 



used approaches and presented a microphysical computational 
model of Knudsen diffusion in random porous media (RPM) 
formed by packed monodisperse spheres. The main transport 
characteristics, such as the mean free path distribution and the 
relative permeability of the porous slab, were calculated for 
a high-porosity medium (porosity was above 65 %) and com- 
pared with that obtained by classical capillary models. Both 
approaches (RPM and capillary models) predict that, in the 
stationary stage of diffusion, the relative permeability of the 
porous layer is inversely proportional to its thickness. This rela- 
tion was observed also in all our experiments. At the same time, 
a clear difference between theoretical and experimental data ob- 
tained for the thin layer has been detected. For example, we got 
an experimental permeability of a monolayer of dust that is sig- 
nificantly higher than the theoretical predictions. We also note 
that the absolute weakening of the gas flux was smaller in the 
experiments. This lack of correspondence between model and 
experiment may be caused by an ill-determination of the bound- 
ary condition in the theoretical models, when the molecules 
scattered back were excluded from further consideration. Obvi- 
ously, new laboratory and computer experiments are necessary 
to clarify the situation and to validate the new, more sophisti- 
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cated theoretical models of cometary nuclei. 

The second focus of this research is the evaluation of the sub- 
limation coefficient. Although the variability of this character- 



istics is confirmed experimentally (see Sect. 3.3 I, the effect has 
not received the attention it deserves in cometary research. The 



publication by |Kossacki et al. ( 1999 1 remains one of a few at- 
tempts to explore the consequences of this phenomenon for the 
transport processes in cometary nuclei. It was shown, for ex- 
ample, that for small heliocentric distances (where the surface 
temperature is above 190K), the reduction of the sublimation 
coefficient leads to a significant growth of the surface tempera- 
ture (up to 20K) and to a change of the temperature distribution 
near the surface. The heating causes trapped volatiles in the 
ice to sublimate more actively than simple models predict. The 
deviation of the sublimation coefficient from unity and the re- 
sulting temperature increase can cause crystallization of amor- 
phous water ice at larger heliocentric distances and/or for larger 
nucleus depth. The increase in temperature associated with the 
damping of the sublimation coefficient leads to an exponential 
increase of the corresponding saturation pressure. We note that 
the covering porous dust layer also induces an increase of the 
surface temperature of the ice (the so-called "cooking effect"). 
Thus, one can expect that both phenomena explored in this pa- 
per can produce multiple effects leading to an increase of the 
gas pressure beneath the porous dust layer. If we assume that 
the gas pressure is the main driver for the breakup of the dust 
crust, such a pressure increase is a prerequisite for the explana- 
tion of dusty outbursts. 
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